Qubits as devices to detect the third moment of current fluctuations 
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Under appropriate conditions controllable two-level systems can be used to detect the third mo- 
ment of current fluctuations. We derive a Master Equation for a quantum system coupled to a 
bath valid to the third order in the coupling between the system and the environment. In this 
approximation the reduced dynamics of the quantum system depends on the frequency dependent 
third moment. Specializing to the case of a controllable two-level system (a qubit) and in the limit 
in which the splitting between the levels is much smaller than the characteristic frequency of the 
third moment, it is possible to show that the decay of the qubit has additional oscillations whose 
amplitude is directly proportional to the value of the third moment. We discuss an experimental 
setup where this effect can be seen. 

PACS numbers: 



I. INTRODUCTION 

A comprehensive understanding of the transport prop- 
erties of mesoscopic conductors can be achieved with the 
study of both the average current and its fluctuations. 
The investigation of shot noise 1-3 has proven to be a valu- 
able tool to determine properties which are elusive to the 
study of current-voltage characteristics. One of the most 
celebrated examples in this respect is probably the mea- 
surement of the fractional charge by means of the study 
of shot-noise in point contacts in the Fractional Quantum 
Hall regime 4 ' 5 . 

In the case of non-gaussian fluctuations, moments be- 
yond the second one are relevant in characterizing the 
transport. In the last few years numerous theoretical 
studies (see Ref.6) analyzed the properties of higher mo- 
ments and of the full counting statistics 7 . In contrast to 
the large theoretical activity, experiments are very diffi- 
cult and only few appeared so far. The first pioneering 
measurement of the third moment, performed by Rculct 
et aP , has been hindered by environmental effects 9 . The 
first three moments of the current fluctuations in a tun- 
nel barrier were measured very recently by Bomze et al. 10 
confirming the Poisson statistics associated with the dis- 
creteness of the charge. Further experimental indications 
on the non-Gaussian character of noise were obtained by 
Lindcll et al. 11 who observed its effects on a Coulomb 
blockade Josephson junction. 

In parallel with the first experiments, and with the 
hope of finding more effective ways to measure higher 
moments, several theoretical papers appeared suggest- 
ing ways to find signatures of non-gaussian noise in the 
non-equilibrium properties of mesoscopic systems used as 
detectors. The first practical way to probe high moments 
of current was suggested by Lesovik in Ref.12. In Ref.13, 
Aguado and Kouwenhoven considered the possibility to 
use a double quantum dot system as detector of high fre- 
quency noise. More recently, Josephson junctions were 



shown to be able to act as detectors of the third 14 and 
fourth moments of current fluctuations 15 . Their use as 
threshold detectors to measure the full counting statistics 
has been discussed by Tobiska and Nazarov 16 , by Leso- 
vik, Hassler and Blatter 17 and by one of the authors 18 . 

Qubits have been already proven very sensitive spec- 
trometers of noise 19,20 . In this work we want to explore 
further the use of controllable two-level systems as noise 
spectrometers and analyze the possibility to employ them 
for the measurement of the third moment as well. With 
this scope in mind we derive a perturbative equation for 
the dynamics of two-level systems in presence of noise 
up to the third order in the system-noise coupling to see 
if, under some circumstances, we can extract some infor- 
mation on the third moment. In general the third order 
effects are masked by the dominant second order ones, 
since they are a result of a perturbative expansion. There 
are, however, situations in which the second-order correc- 
tion vanishes and therefore the third order is the leading 
contribution. We will show that, in the usual Rotating 
Wave Approximation (RWA) of the system equations of 
motion, the contribution of the third moment is a small 
correction to the dominant effect of the second moment 
(and hence difficult to measure). A treatment beyond 
RWA is therefore needed and it leads to the presence of 
additional effects solely due to the third moment of cur- 
rent fluctuations. 

The paper is organized as follows. In the next Sec- 
tion we introduce the model. We then proceed to derive 
the Master Equation for the reduced dynamics of the 
quantum system up to third order in the coupling with 
the environment. The reduced dynamics of the quantum 
system will also depend on the third-order correlations 
of noise and therefore it may act as a detector of these 
higher order moments. In Section III A we concentrate 
on the case in which the quantum system is a two-level 
system and show that the presence of the third order 
may induce coherent oscillations in the ground state pop- 
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ulation of the quantum system. Furthermore, we show 
that the amplitude of these oscillations is proportional to 
the three-point correlator of the fluctuations. In Section 
IIIB we discuss the case where the two- level system is 
subject to an external microwave field. In particular we 
discuss how Rabi oscillations can be influenced by the 
presence of the third order noise. The motivation here is 
to lower the frequency of the coherent oscillations into a 
regime which would be more accessible to experiments. 
Possible experimental setups where these effects can be 
measured are discussed in Section IV. In the same Sec- 
tion we analyze various complications which may emerge 
in the actual measurement. In particular we consider the 
case of a DC-SQUID as a third-moment detector. Sec- 
tion V is devoted to a summary of the results and pos- 
sible perspectives of this approach in measuring higher 
order current fluctuations. Recently a similar detection 
scheme has been discussed in Rcf.2I; our approach is dif- 
ferent in spirit and we will point out the difference with 
Rcf. 21 where the Rotating Wave approximation is taken 
for granted 22 . 



II. DYNAMICS OF THE QUBIT 

A. The model 

In this first section we discuss the general formalism 
which will be used in the remainder of the paper. The 
setup under consideration is composed of a quantum sys- 
tem S weakly interacting with a quantum bath B. As 
explained in the Introduction the quantum system will be 
used to investigate the properties of the bath which for 
example may be another nanostructure (a tunnel barrier, 
a point contact,...) biased at a fixed voltage and of which 
we want to study current fluctuations. The Hamiltonian 
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of the total system S + B can be written as follows: 

H T = H S + H B + V, (1) 

where H$ and Hb are respectively the free Hamiltonian 
of the system and of the bath. The interaction potential, 
V, is chosen to be of the form 



V = g^N a ®Qc 



(2) 



In the definition of V, g is an adimensional coupling con- 
stant and N a and Q a are operators of the bath and sys- 
tem respectively. The interaction is chosen to be weak so 
that the dynamics of the reduced density matrix of the 
system can be obtained by a perturbation expansion in 
g. The procedure is well known and described in various 
textbooks 24 . It is typically performed up to second order 
in the coupling g; here we do a step forward and go to 
the next order in the coupling. As we focus our atten- 
tion on the study of the time evolution of the system S 
in presence of a stationary bath, we have [pb,Hb] = 0, 
ps being the bath density matrix. Moreover we assume 
the dynamics of the whole system to be Markovian. This 
means that at each order of perturbation theory we can 
neglect all the terms that are non-local in time. 



B. Third order Master Equation 

The time evolution of the reduced density matrix of the 
system in the interaction representation is described by 
the following third order equation (the steps leading to 
the master equation are standard 24 and we do not repeat 
them): 



pi = Tn 



■ j£ jTdf [V(t), [V(t% pj(t)p B ] ] + ^ jHdf J* dt" [V(t), [V(t>), [V(t"), pi(t)p B }]} I , (3) 
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where we denoted respectively with pi and V(t), the den- 
sity matrix of the system and the interaction potential 
in the interaction representation. In deriving Eq.(3) we 
made the further assumption that (V) = Ttb\pbV\ = 0, 
where Tr B denotes the trace over the bath degrees of 
freedom. Taking the matrix elements of (3) between two 
eigenstates of the Hamiltonian of the system, after some 
algebra, we obtain the following third order master equa- 
tion for the density matrix of S: 



Pi, ran = ^ P I ,lk,D Unml^' 



(E m -E n -E k +Ei)t 



(4) 
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where we set H$ \m) — E m |m), pj, mn = (m\pi\n) and 
the third order relaxation matrix, Dknmi , is given by the 
sum of two contributions: 



Dknml — 2 ^kn 



_ l JLr 

n 



(5) 



In the previous equation Rknmi is the second order relax- 
ation matrix and Cknmi is a third order correction crucial 
to our treatment. Here we omit to show the explicit ex- 
pression of Rknmi as it is well known and it can be found 
in textbooks 24 . 

The third order kernel Cknmi can be written as follows: 



Cknmi — Aknml ^-Irank 
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mlkn 
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nklm J 



(6) 
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with 



*-knml 



EE 

a/37 i 



E QmiQijQ]l S knF a /3j(uJli,C0ij) ~ QmiQllQkn F aP-f(ui m , U>mi) (7) 
3 

Bmlkn = E E [QtlQLQfnGap-fiuJlm + UJ ik ,U m l) - QZiiQuQln G ^{^U + ^nk,^l) (8) 
a/37 » 

and where 7iw mn = £" m — E n . The functions F Q/ 3 7 , G Q( g 7 are the three-point correlators of the noise operators: 

roc rt-t 

FafriwuUi) = dh dt 2 (N a {t 1 )N l3 (t 2 )N 1 (Q)) B e l ^e 1 ^ (9) 
Jo Jo 

/>oo />t± 

G a0J (u; 1 ,W2) = dh dt 2 (N J (0)N a (t 1 )N^t 2 )) B e^e^ . (10) 
Jo Jo 



We used Q„„ instead of (m Q a nj, to denote the 
matrix elements of a system operator in the Schrodingcr 
picture. 

The average = Ttb[pb ■■■] is taken over the density 

matrix of the bath. Note that we do not need the bath to 
be in equilibrium but we do assume that it is stationary. 
An example is the noise generated by a non-equilibrium 
current in a voltage biased tunnel junction. 

Equation (4) is quite general and many specific cases 
can be studied starting from it. In the Rotating Wave 
Approximation, i.e. neglecting oscillating terms in the 
sum on the right-hand side of Eq. (4) , one can recover the 
result of Ref. 21. In this case the presence of the third 
order causes simply a small correction of the second order 
transition amplitudes and therefore it might be difficult 
to detect in the presence of a large background due to the 
second order contribution. In the following sections, we 
will analyze in more detail some special cases in which 
the effects of the third order relaxation matrix can be 
well characterized and distinguished from the second or- 
der. In particular, we will study how the third order 
contribution affects the decay and the Rabi oscillations 
of a two-level probe quantum system. 



III. RESULTS 

We now specialize to the case in which the probe is a 
two-level quantum system. We assume that the effective 



r 



Hamiltonian of the system in the presence of noise has 
the form 



H e ff = - fruod-z + N T (t)a x 



(11) 



when expressed in the eigenbasis of S, iVy(t) is the noise 
operator and &i are the Pauli matrices. Moreover we 
make the hypothesis that the relevant frequencies of 
the noise source are much larger than the level split- 
ting ujq. We thus neglect the frequency dependence of 
the third order correlators on scale up to luq. Conse- 
quently in the calculation of the third order coefficients 
of the relaxation matrix, we set F(u>i,uj 2 ) ~ F(0,0) and 
G(ui,u 2 ) ~ G(0,0), if u)\,uj 2 ~ ivq. In the following we 
will comment on these assumptions. 



A. Relaxation in presence of non-gaussian noise 



In case of a two-level system the third order master 
equation, Eq.(4), in the Schrbdinger representation, re- 
duces to 



J 



pn = (Dun - D 211 2) pn + (D 1112 + D 2111 )Re[p 12 ] - i(D 1112 - D 2111 )lm[p 12 ] + D 2112 (12) 
P12 = (-D1211 - -D2212) pu +-D2212 - i OD1212 - Ami - ico )lm[p 12 ] + (-D2211 + D 1212 + iu; )Re[p 12 }. (13) 

I 



where Tiloo = E 2 — E\. The different elements of the third the definition given in the previous paragraph, Eqs. (5)- 
ordcr relaxation matrix, Dknmh carL be calculated using 



(10). Within our hypothesis, the only nonzero second 
order contributions are 



P,,W 



-D2112 
£>im 
-D1212 
-D2211 



2 roo 

= 9 dt > (N T (t')N T (0)) e^*' = W 21 (U) 

*l J —OO 

= -D 2 ii2(-wo) = -W 12 (15) 

„2 /-oo 

- fa / dt'<{JVT(f),JV r (0)})e-«'*' (16) 

Al Jo 

(17) 



-D 



1212- 



where we have introduced the second order transition 
rates, W\ 2 and W 2 i- 

Note that, due to our transverse coupling assumption, the 
third order contribution to the previous matrix elements, 
Eqs. (14)-(17), is zero. The other coefficients of the re- 
laxation matrix are of the third order in the coupling 
constant, g. In the limit of a flat spectrum all these 
elements can be defined as follows, using only one inde- 
pendent parameter 



12 = 



D 



2212 



^A (3) D 2111 = D 1211 =iA (3) (18) 



The third order coefficient, A^ 3 \ is real and it can be 
written a sum of time-ordered products as follows 




AO) = S 

h 



3 , dti I dt 2 



{N T (h)T (N T (t 2 )N T (0))) + 



+ -{T(NT(ti)N T (t 2 )NT(0)) 



(19) 



where T and T denote respectively the time-ordering 
and the anti-time-ordering operator. 

Third moment fluctuations can be measured by mea- 
suring the probability that the system is in the ground 
state once it was initially prepared in the state 

|*(* = 0)) = -^(|1) + |2)). 

Indeed the ground state population as a function of time 
can be easily calculated from the integration of Eqs. (12)- 
(13) 



Pn(t) = 



W21 



W 12 + W 2 i 



+ e -^ +w ^ t [A + 2B cos{uj R t)] . 

(20) 

In the previous equation, we have introduced the renor- 
malized frequency 

u>% = wg + uj a lm [D 1212 ] - i(Wi2 + Vt^2i) 2 (21) 

and the coefficients A and B which are defined by the 
following equations 



1 14^12 - W 21 
A = - - " 



A (3) 



B = 



2 W 12 + W 21 2uj r 
4^ + °^- 



+ 0( 5 4 ) 



(22) 



(23) 



Figure 1: Ground state population as a function of time and 
its Fourier transform calculated for the following values of 
the parameters: W12 + W21 = uog 2 , W12 — W21 = 0.05 uiog 2 , 
A (3) ~ 0.3u;o£ 3 , Wfl = u> , g 2 = 0.1 . 



The presence of the third order, or more precisely, the 
presence of nonzero odd moments of noise fluctuations, 
induces measurable effects in the dynamics of the probe 
quantum system. As one can see from Eqs. (20)- (23), it 
induces coherent oscillations in the ground state popula- 
tion of amplitude proportional to the third order param- 
eter, A (3) . 

In Fig.l we show the ground state population as a 
function of time and its Fourier transform. As one can 
see pn(t) is given by the superposition of two terms: 
a damped exponential whose asymptotic value is fixed 
by the ratio between the two transition amplitudes, W\ 2 
and W21, and a damped cosine term proportional to the 
third moment. The structure of pn(t) can be analyzed 
by studying its Fourier transform (shown in Fig.l lower 
panel) defined as 



/» -1 r 

j — c 



dte^ [/(*)-/(* = 00)] 



The zero frequency peak is related to the second order 
non-oscillating contribution, while the smaller peak at 
frequency ujo is a pure third order effect. In absence of 
third moment the time dependence of the ground state 
population would be simply described by a damped ex- 
ponential and no third order peak would appear in the 
Fourier transform at lo = loq. 
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The assumption that the noise couples to the sys- 
tem only through a x (transverse coupling) was crucial 
in our analysis to separate the second and the third or- 
der contributions in different elements of the relaxation 
matrix. This assumption can be relaxed, by introduc- 
ing in the Hamiltonian a longitudinal term of the form: 
Vt(t) = gLNL(t)a z , provided that the two noise opera- 
tors Nx(t) and N^it) can be considered as uncorrelated 
and that Vl(£) is weak. In this case the final result is 
essentially the same except for a redefinition of order g\ 
of the transition amplitudes and of the renormalized fre- 
quency UJ R . 

B. Effects of a microwave field 

As shown in the previous section the presence of odd 
moments in the current fluctuations has a distinct signa- 
ture in the oscillations of the ground state population in 
the case of transverse coupling to the noise, Eq.(ll). 

However, the actual measurement of these oscillations 
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Due to the assumption of transverse coupling to the 
noise source and of the frequency independence of the 
third order correlators, the different coefficients of the 
third order relaxation matrix are given by Eqs.(14)-(19). 
Note that setting the coefficients Dknmi and the longi- 
tudinal microwave contribution Ml to zero one easily 
recovers Rabi theory. In this case, solving the eigen- 
value equ ation, one finds t he known Rabi frequency: 
uiRabi = ^M%/4+(n-uj) 2 . 

We first discuss the outcomes of a numerical integra- 
tion of Eqs.(25)-(26). The coupling constant and the 
renormalized frequency are the same in all the figures: 
u R = ujq , g 2 = 0.1. 



In Fig. 2 we show pn{t) and its Fourier transform in pres- 
ence of a weak transverse microwave field. The structure 
superimposed to the damped Rabi oscillation can be 
better understood by looking at the Fourier transform. 
We see indeed four peaks. The zero frequency peak 
corresponds to a pure damping term. At the detuning 



can be very difficult as their characteristic frequency, 
lor ~ ujq, is typically of the order of 10 GHz and the 
time resolution required to follow such oscillations in de- 
tail is hardly accessible. In this Section we discuss a 
generalization of the the case discussed before to account 
for the presence of an external microwave field. Our aim 
is to clarify under which conditions a microwave field can 
shift the third order peak to a lower value fixed by the 
detuning frequency. 

In presence of microwaves, the dynamics of the quan- 
tum system can be described by the effective hamilto- 
nian: H e ff = Hs + V + M. The effect of the applied 
field leads to the term M = O cos(fii) where O is a sys- 
tem operator which, quite generally, can be expressed in 
the form 

0=| (M L a z + M T a x ) . (24) 

The corresponding Master Equation for a two-level 
quantum system in the presence of a microwave field is 
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frequency, — u)q, we see a large Rabi peak whose am- 
plitude does not essentially depend on A' 3 ) and whose 
width is fixed by the relaxation rate T = W21 + W\2- 
At frequency luq we find the contribution arising from 
the third moment fluctuations which is essentially not 
modified by the presence of a weak transverse microwave 
and that depends on A^ 3 ) linearly. The last small peak 
at il + ujq does not originate from the third moment but 
is due to non-secular terms already present in second 
order in the coupling to the environment. 

In Fig. 3 we show our results in case of a strong longitu- 
dinal field. In order to suppress the second order effects, 
we set W12 = W21; then, in the absence of third-order ef- 
fects, one would simply have a constant ground state pop- 
ulation equal to its asymptotic value pn(t) = = 0.5. 
The non-trivial time-dependence in the ground state pop- 
ulation is therefore completely related to the presence of 
the third moment fluctuations. In the case of longitudi- 
nal field in Fig. 3 (lower panel), there are two peaks of 
non-negligible amplitude in the Fourier transform at the 
frequency luq and at the detuning frequency, respectively. 
The u>o peak is the peak present also in the absence of 
the microwaves (both its amplitude and position are es- 
sentially not affected by the presence of the microwave 



(Dim - D 2 n2)pii + (D1112 + D 2 iii)Re[pi 2 ] + D 2 n 2 - i[D 1U 2 - D 2m - iM T cos(Qi)]lm[pi 2 ] 

(25) 

[Dim - D2212 + iM T cos(ftf)] pi\ + [M L cos(Oi) - ui - i (D m2 - D 22 n)] Im[pi 2 ] 

+ [D2211 + Di2i 2 - i (M L cos(m) - lo )] Re[p 12 ] + D 22 i 2 - iM T cos(m). (26) 




p„(t) 




Figure 2: Ground state population as a function of time (up- 
per panel) and its Fourier transform (lower panel) calculated 
in presence of a weak transverse microwave field. We chose 
W12 = W21 = 0.5g 2 uj , A (3) = 0.6pV) ~ 0.002, Q = 1.3w , 
M T = O.Icjq, and M L = 0. 



Figure 3: Ground state population as a function of time (up- 
per panel) and its Fourier transform (lower panel) calculated 
in presence of a strong longitudinal microwave field. The pa- 
rameters are W12 = W21 = 0.5 g 2 io , A (3) = 0.3 g 3 co ~ 0.001, 
n = 1.3w , M L = 0.8w , and M T = 0. 



field). The second peak, located at frequency Q — u>q, is a 
combined effect of the third moment fluctuations and of 
the microwave field; its amplitude is directly proportional 
both to the value of A^ 3 ) and of M L . The position of this 
peak is determined solely by the detuning frequency and 
it is not affected by the amplitude of the microwave field, 
M L . 

In the case of pure longitudinal field an approximate 
analytical solution of Eqs.(25)-(26) can be found. Up to 
third order in the coupling constant, g, we obtain 

pn(t) - a + be-^ +w ^ t + 

00 (M \ 
+c ^ Jk{-^-)sm((uJo-kQ)t + 6k) ■ 

(27) 

In the previous equation Jk(z) is the k-th Bessel func- 
tion 25 , the phases 9 k and the real constants a, b and c 
are defined as 



0k 



arctan 



12 



W 21 



w 21 



f W21 

AW 
40 



OO 



(28) 

J fe (^)sin(0 fe )(29) 



W 21 



2 W 12 + W21 

A(3) 

4fl 



(30) 
(31) 



The peak at a>o is associated to the k — contribution of 
the sum while the detuning frequency peak and the peak 
at frequency Q + lj are related respectively to the k = 
±1 -terms. The small contribution connected with the 
k = 2 Bessel function is also visible in Fig. 3 at frequency 
uo = 2Q — uoq. 

In the hope to use the method discussed in this pa- 
per for the diagnostics of the third-moment of current 
fluctuations it is useful to analyze the amplitudes of the 
different peaks in some detail. In Fig. 4 the height of the 
peak at loq is shown as a function of A^ 3 ) . In Fig. 4a we 
show the results in case of weak fields. In this case, the 
presence of the longitudinal or the transverse field does 
not affect the height and the position of the third order 
peak. In Figure 4b we display the results in case of strong 
fields. As one could expect from Eqs.(25)-(26), a strong 
transverse field masks completely the third order depen- 
dence of the Wo peak; on the other hand, even a strong 
longitudinal microwave field does not essentially modify 
the height and the position of the third order peak at 
frequency u>o- Note that the range [0, 0.03] for A^ 3 ) /u>o is 
chosen so that the ratio between second and third order 
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0.005 0.01 0.015 0.02 0.025 co/co, 
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0.005 0.01 0.015 0.02 0.025 co/co, 
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0.005 0.01 0.015 0.02 0.025 A <3) /(|) 



Figure 4: Amplitude of the peak at u>o as a function of 
the third order correlator A^ 3 ' in absence of microwave field, 
(solid line), and in presence of transverse (dotted line) and 
longitudinal (dot-dashed line) fields, a) Weak microwave 
fields respectively with Mt = 0.1 wo or Ml = O.lwo. b) 
Strong microwave fields respectively with Mt = 0.8 u>o or 
Ml = 0.8^0- The other parameters are the same in both the 
cases: W12 + W21 = g 2 coo, W12 — W21 = 0.05 coog 2 , = 1.3a>o- 



contribution varies between and g. 
Figures 5 and 6 are devoted to the study of the amplitude 
of the peak at the detuning frequency, f2 — u>o, in case of 
pure longitudinal field. In order to compare the numer- 
ical results shown in these figures with the approximate 
analytical solution (27), we now give the explicit expres- 
sion of the k = 1 term, which is responsible for the peak 
at the detuning frequency. This term can be rewritten as 



1)^ 



cJl (^) sin (A*- 

A< 3 >M £ (VFi 2 + W 21 ) sin (At) - A cos(At) 
8fl (W12 + W21) 2 + A 2 1 J 



where we set A = il — ui n and in the last step we kept 
only the linear term in the field amplitude. In Fig. 5(a) 
and Fig. 5(b) we display the amplitude of the peak, re- 
spectively, as a function of A^ 3 ) and of M L . As one could 
expect, the amplitude of the peak is proportional to A' 3 ); 
moreover, as one can see in Fig. 5b, the linear approxima- 
tion is fulfilled also in case of strong fields. In Fig. 6, we 
show the amplitude of the peak as a function of the de- 
tuning frequency; again the result is as expected based 
on Eq.(32). 



Figure 5: Amplitude of the peak at the detuning frequency 
Q — uio as a function of of the amplitude of the microwave field 
Ml (a) and of the third order correlator A' 3 ^ (b). Values of 
the parameters: W12 = W21 = 0.5uiog 2 , = 1.3 wo, Mr = 0. 
In (a) M L = 0.8 luq while in (b) A (3) = 0.01 u) . 




Figure 6: Amplitude of the peak at the detuning frequency 
as a function of the — cjo under longitudinal microwave 
coupling. Values of the parameters: W12 = W21 = g 2 loo, 
M T = 0, M L = 0.8u , A (3) = 0.01 w . The inset is a magni- 
fication in the region Q ~ loq. 



IV. EXPERIMENTAL CONSIDERATIONS 

Starting from the results presented in the previous sec- 
tions, we would like to discuss an experimental protocol 
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for the measurement of third order noise using a two-level 
probe quantum system. The experimental realization of 
this kind of measurements is rather delicate. The ba- 
sic reasons for this is that a set of inequalities has to be 
satisfied. First, the level spacing has to exceed the tem- 
perature in the experiment, T, in order to avoid ther- 
mal excitations: Hujq 3> ksT. Second, the time reso- 
lution in the experiment, St, has to be good enough to 
follow the detuned coherent oscillations at angular fre- 
quency f2 — u>o, i.e., St <C |n — ^ol _1 - Yet the oscillation 
in the ground state population has to be measurable, 
which means that |(Q — u)o)/wq\ should not be too close 
to zero. This condition is determined by the resolution 
in measuring the population variations: for very small 
|(0 — u>o) /ljq |) -Po(0 — wo) is significantly suppressed as 
demonstrated in Fig. 6. Collecting these conditions we 
have 



noise 
source 



\(n-L) )/u \ < 



i 



u>o5t ksTSt 



(33) 



As to a concrete realization, one may employ a hysteretic 
DC-SQUID in the tunnelling regime 23 . The strength of 
the method lies in the high contrast in resolving level oc- 
cupations: tunnelling rate from the excited state is typi- 
cally two to three orders higher from the excited state as 
compared to that from the ground state. Measurement 
of this decay is straightforward by observing the switch- 
ing statistics, i.e., the measurement is typically repetitive 
in nature. Occupation probabilities of order 0.1 or even 
below are measurable with adiabatic detection pulses of 
St ~ 1 ns duration. Measurements are typically car- 
ried out at T ~ 30 mK. Taking this temperature, the 
condition at the right end of Eq. (33) then states that 
ujo ^S> 3 • 10 9 s _1 . Typical level separations (plasma fre- 
quencies) of Josephson junctions are in the range 1 GHz 
< ojq/2tt < 100 GHz; thus these values are compatible 
with the operation temperature. With ujq ~ 10 10 s _1 and 
St = 1 ns, we match the frequency vs. temperature con- 
dition with a margin of factor three. The other critical 
condition in Eq. (33), |(0— u)q)/u>q\ <C (uoii) -1 , can then 
be matched by requesting |(Q — Wq)/ujq\ <C 0.1. Since the 
maximum of Pq is obtained at \Q/u>o — 1| ~ 0.05, we no- 
tice that both the inequalities can be satisfied, although 
barely. 

The remaining questions then concern the coupling of 
the noise to the detector. As we have already pointed out, 
in order to see oscillations in the occupation probability 
due to purely third order effects, one needs to couple the 
noise source to a x and the microwaves to a z . In case of a 
DC-SQUID detector, this means that one should couple 
the noise source through the current and the external mi- 
crowave field through the flux. A schematic diagram of a 
possible measuring apparatus is shown in Fig. 7. The de- 
tector is constituted by a DC-SQUID of negligible induc- 
tance formed with two identical Josephson junctions of 
critical current Iq and of capacitance Co/2 biased by ex- 
ternal flux, and current, /. The non-equilibrium noise 
source, which can be a tunnel junction or another nanos- 



SQUID 
detector 




Figure 7: Schematic diagram of a possible experimental setup 
to detect third order fluctuations 



tructure, induces time dependent fluctuations in the bias- 
ing current, SI(t). Finally, the external microwave field 
is coupled inductively to the SQUID ring and leads to 
monochromatic flux fluctuations, S$(t). The effective 
two level hamiltonian of this system including flux and 
current fluctuations can be written as follows 



H 



TlLOr, 



Tiuj p [N T (t)a x - N D (t)a z 



(34) 



The operators iVr(i) and Noit) are defined by the fol- 
lowing equations 

N T {t) = t I —^ + (t I -2\d$)—^- (35) 
N D {t) = (d^ + 2\t$)^P^ + 2\t I ^P- (36) 



Io 



where <&o is the elementary flux quantum. In the previous 
equations we denoted with A and aj p , respectively, the 
SQUID anharmonicity and the SQUID plasma frequency: 



W7, 



2tt 



Vlf^P (37) 



with I c = 2I cos || and r) = h* (C / c )" 
Moreover we have introduced the adimensional parame- 
ters ti, t<£, 

1 I ( I 2 \~* I 7T$ , 7T$ 

t/ = ^l 1_ ^J ' * $ = 7o* /tan V d#=tan «S- 

(38) 

As it is clear from Eqs.(34)-(38) a transverse coupling 
to current fluctuations and a longitudinal coupling to flux 
fluctuations (i.e. to the microwave field), can be simul- 
taneously realized only in case of zero or very small DC 
component of the external biasing current. Experimen- 
tally this condition could be obtained subtracting the DC 
component of I(t) by means of a superconducting line, as 
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it was shown in Ref.26. In this case the SQUID potential 
becomes harmonic and one should also take into account 
the possibility to have transitions to higher levels. Any- 
way, if initially only the first two levels are occupied, this 
effect is of order higher than the third in the system- noise 
coupling g. 

Assuming that the transverse coupling condition is ful- 
filled, that is / << I c , we now evaluate the amplitude 
of the third order oscillations for the detector's scheme 
presented in Fig. 7. 

Comparing Eq.(34)-(38) to the model hamiltonian in 
Eq.(ll) and using the definition given in equation (19), 
we rewrite A 1 - 3 ) as follows 



A (3) = [ 

1 



dti / dt 2 



(6I(h)T(SI(t 2 )6I(0))) + 



f(SI(h)SI(t 2 )SI(0)) 



(39) 



Using the results derived by Salo, Hekking and Pekola 
in Rcf.27 within the framework of scattering theory 2 , we 
can easily obtain an explicit expression of A^ 3 ) in terms 
of the transmission eigenvalues, T n , and of the voltage 
bias, V, across the junction. 

In particular, in case of energy independent scattering in 
the limit of zero temperature of the noise source, we have 



A< 3 ) = 



3 \rjlc 



eU ^ 



T„(l-T n )(l-2T„); (40) 



as one can see, in this limit, A' 3 ) is proportional to the 
usual third cumulant of current statistics 7 ' 28 . Eventually, 
using Eq.(40), we can also check the validity of the per- 
turbative hypothesis and give a rough estimate of the ra- 
tio between the third and the second order contributions 
to the qubit dynamics. In the limit of zero frequency and 
zero temperature we have: 



A( 3 ) 2_F 3 



3 9 Y 2 



(41) 



CONCLUSIONS 



In this paper we analyzed the possibility to use solid 
state qubits as detectors for higher moments of current 
fluctuations. We showed that in some cases there are dis- 
tinct features, due to the non-secular terms in the Master 
equation, solely related to the presence of the third mo- 
ment of current fluctuations. This may be a very inter- 
esting circumstance as usually these additional effects are 
masked by the large background coming from the noise 
(second order cumulant in the fluctuations). After hav- 
ing derived the general form of the Master equation up to 
the third order in the coupling between the environment 
and the bath we considered in some detail a two-level 
system coupled to a noise source. Indeed we found that, 
in the presence of purely transverse noise, the population 
in the ground state oscillates at a frequency wo if the two- 
level system is initially prepared in a superposition. The 
difficulty of measuring these high frequency oscillations 
can be alleviated by applying a microwave field. In this 
case the oscillations associated to the third-moment are 
pushed down to the detuning frequency Q — ujo ■ 

A possible experimental implementation of this scheme 
of detection has been discussed in Section IV. As a two- 
level system (the detector) we considered a DC-SQUID 
and discussed the range of applicability of the scheme. 
Combining the narrow margins in experimental parame- 
ters and the rather unfavourable coupling of noise to the 
detector, it is obvious that measurement of the effects 
predicted here is not straightforward using a DC-SQUID 
as a sensor. It remains to be analyzed if other control- 
lable two-level systems (charge qubits for example) may 
be more suited as detectors. Nevertheless we find inter- 
esting the existence of features entirely due to the higher 
moments of current fluctuations. 



where g = eiv p /(r]I c ) ~ (Co J c )~* ■ (2.1 1CT 6 (A • farad)*) 
is an adimensional coupling constant and F 2 and F 3 are 
the Fano factors of the second and of the third order: 
F 2 = £„Tjl-T^/£ n T rl and^3 = E„7Ul-T„)(l- 
2T , „)/]^; rl T , „. In deriving Eq. (41) the well-known rela- 
tion between the second order transition amplitude and 
the the Fano factor F 2 has been used, see, for example, 
Refs.1,13. 
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